Disentangling the electronic and phononic glue in a high-Tc superconductor 
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pH ' Unveiling the nature of the bosonic excitations that mediate the formation of Cooper 

^ I pairs is a key issue for understanding unconventional superconductivity. A fundamen- 

fj , tal step toward this goal would be to identify the relative weight of the electronic 

I ' and phononic contributions to the overall frequency (SI) dependent bosonic function, 

i^ ' n(Sl). We perform optical spectroscopy on Bi2Sr2Cao.92Yo.o8Cu208+i crystals with si- 

r-{. multaneous time- and frequency-resolution; this technique allows us to disentangle 

jy^ ' the electronic and phononic contributions by their different temporal evolution. The 

strength of the interaction (A~l.l) with the electronic excitations and their spectral 
C^ , distribution fully account for the high critical temperature of the superconducting 

phase transition. 

i-rt ' 

r^ \ Lattice vibrations [11] and excitations of electronic origin, like spin or electric polarizability fluctuations [3] and loop 

Q ■ currents'3], are generally considered potential mediators of Cooper-pairing in the copper-oxide high-temperature 

O ] superconductors (cuprates). The generic interaction of fermionic quasiparticlcs (QPs) with bosonic excitations 

is accounted for by the bosonic function 11(17) (usually indicated as a'^F{n) for phonons and Px{^) for spin 

fluctuations), a dimensionless function that depends on the density of states of the excitations and the strength 

J> , of their coupling to QPs. Because both the energy dispersion and lifetime of QPs are strongly affected by the 

OO ' interactions, signatures of QP-boson coupling have been observed in experiments that probe the electronic properties 

OC \ at equilibrium. The ubiquitous kinks in the QP dispersion at ^70 meV, measured by angle-resolved photoemission 

^^ • spectroscopy (ARPES)[1], have been interpreted in terms of couphng to either optical Cu-0 lattice modes|3,@| or spin 

^^ \ excitations [7||. Inelastic neutron and X-ray scattering experiments found evidence for both QP-phonon anomalies !8| 

CO , and bosonic excitations attributed to spin fluctuations [7|, |9| and loop currentsllOJ. Dip features in tunnelling 

^^ ' experiments have been used to alternatively support the scenarios of dominant electron-phonon interactions |ll| 

2_| . or antiferromagnetic spin fluctuations fl2|. The frequency-dependent dissipation of the Drude optical conductivity, 

. . ' ct(w), measured by equilibrium optical spectroscopies, has been interpreted [l3l4l5| as the coupling of electrons to 

^ , bosonic excitations, in which the separation of the phononic and electronic contributions is impeded by their partial 

k>( ■ coexistence on the same energy scale (<90 meV). 

C^ ■ We disentangle the electronic and phononic contributions to n(J7) through a non-equilibrium optical spectroscopy, 

in which the femtosecond time-resolution is combined with an energy-resolution smaller than 10 meV, over a 
wide photon energy range (0.5-2 eV). Our approach is based on the widely-used assumption p^ \V(\ that, after 
the interaction between a superconductor and a short laser pulse (1.55 eV photon energy), the effective electronic 
temperature (Te) relaxes toward its equilibrium value through energy exchange with the different degrees of 
freedom that linearly contribute to 11(17). In a more formal description, the total bosonic function is given by 
n(ri)=nbo(17)-|-nscp(^)+niat(f7) where Ilbc refers to the bosonic excitations of electronic origin at the effective 
temperature Tbe, IIscp to the small fraction of strongly-coupled phonons (SCPs) at Tscpflal and Iliat to all other 
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FIG. 1: Non-equilibrium reflectivity. The relative reflectivity variation SR/R{uj)= [R{uj,T + 5T)-R{uj,T)]/ R, calculated for 
the incremental variations STe and 5Ti, in Eq. 1 and using the n(aj) obtained from the fit to the equilibrium measurements, 
is reported in the Fig. SI. The change of sign of SR{ijj) /R in the quasi-thermal scenario (grey line, (5Te=5Ti, = l K) and the 
maximum 5R{u}) /R in the non-thermal scenario (red line, JTe=5 K, (5Ti,=0) coincides with the dressed plasma frequency, 
f2p~l eV. The spectral regions probed by the Optical Parametric Amplifier (OPA) and supercontinuum techniques [Ig] are 
marked by the bars with different colors. 



lattice vibrations at Tiat. As the term that couples the rate equations[lq| for To and T\, (with b=be,SCP,lat) is 
G'(nb,Tb,re)/Cb [G being the functional described in flS] and Cb the specific heat of the bosonic population 6), each 
subset of the bosonic excitations is characterized by different relaxation dynamics on the femtosecond timescale. The 
most convenient systems for such an experiment are the hole-doped cuprates close to the optimal dopant concentra- 
tion needed to attain the maximum critical temperature Tc, in which the total 11(0) is maximum [15l]. Despite the 
magnitude of n(r2), vertex corrections beyond Eliashberg theory can be reliably neglected [9|,[l5[ in this doping regime. 

The total bosonic function 11(17) is directly determined by fitting an extended Drude model[l^ and a sum of 
Lorentz oscillators accounting for the interband optical transitions in the visible region[l9J to the equilibrium 
dielectric function of optimally doped Bi2Sr2Cao.92Yo.o8Cu208-i-(5 (Y-Bi2212) high-quality crvstals |20| (T 'c=96 K), 
measured at T=300 K by conventional spectroscopic elliDSometrv|2l| . Assuming a histogram-like form[l^ (see Fig. 
SI) and imposing an upper limit of 1 eV, the extracted n(i7) is characterized by: i) a low-energy part (up to 40 meV) 
compatible with the coupling to acoustic phonons[22] and Raman-active optical phonons involving c-axis motion of 
the Cu ions [231; ii) a narrow, intense peak centered at ^60 me V, attributed to the anisotropic coupling to either 
out-of-plane buckling and in-plane breathing Cu-0 optical modesp] or bosonic excitations of electronic origin such 
as spin fluctuations [7]; iii) a broad continuum extending up to 350 nieV il3nl5i] , i.e., well above the characteristic 
phonon cut-off frequency (^-^90 meV) . 

The key point to extend this analysis to non-equilibrium experiments is that the electron self-energy, I](aj,i), 
entering in the calculation of ^(w, i)[18], can be factorized into[2J]: 



s(^,r) 



n{n)L{uj,n,T)dn 



(1) 



where L {uj,Q,T) is a material-independent kernel function accounting for the thermal activation of the bosonic 
excitations and of the QPs. The kernel function L (a;,17,Te,b)= -2m[N {n,Tb) +l/2]-h*(l/2-hi(0-0')/27rre)-*(l/2- 
i{n+il,')/2TrTc), where '^ is the Digamma function obtained by integrating the Fermi-Dirac functions and N{i^,T}^) 
the Bose distribution at temperature Tb, can be decomposed into different terms depending on the electronic (To) 
and bosonic (Tb) temperatures. The independent variation of To.b is expected to induce different modifications of the 
dielectric function. Figure 1 shows the expected relative variation of the reflectivity, i.e., SR/R {uj, T)= [R (oj, T + 5T)- 
R(u!,T)]/R{uj,T), in the quasi-thermal {STc—ST\,>0) and non-thermal {STc>0, ST\,=0) scenarios. Phenomenologi- 
cally, in the first case the reflectivity variation is dominated by the increase of the QP-boson scattering, corresponding 
to a broadening of the Drude peak, while in the second case the decoupling between the QP and bosonic distributions 
can be rationalized in terms of a small increase of the plasma frequency without any change in the scattering rate. 
The difference between the two cases is more significant in the spectral region close to the dressed plasma frequency, 
Op~l eV, i.e., an energy scale much higher than the energy scale of the bosonic function. 

Time-resolved reflectivity measurements in the 0.5-2 eV photon-energy range '2511 have been performed at T=300 
K on the same crystals|2(| (OP96) used for the equilibrium optical spectroscopy. The SR/R{u},t) two-dimensional 
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FIG. 2: Non-equilibrium optical spectroscopy. Time- and frequency- resolved reflectivity data carried out on the OP96 sample 
at T=300 K in the OPA (A) and supercontinuum (B) probe optical regionlla]. The white curves represent the time-traces at 
1.55 eV and 0.65 eV photon energy. (C) Energy-resolved spectra for OP96, measured at fixed time i=100 fs. The black lines 
are the maximum and minimum 5_R/_R(aj,i=100fs), determined through the fitting procedure described in the text and in [3| 
and acounting for the experimental uncertainty in the laser fluence and spot dimensions. Most of the values of (5_R/_R(i=100fs) 
measured on OP96 at each wavelength, fall between the two lines. (D) Time-resolved trace for OP96 at fixed 1.55 eV photon 
energy. The black line is the best fit to the data obtained though the fitting procedure described in jig ]. 



matrix is reported in Fig. 2A and B, along with the time-traces at ^1.5 and 0.7 eV photon energies (white curves). 
After the pump excitation at i=0, the temporal dynamics above and below ilp are very similar, exhibiting a relaxation 
dynamics of about 200 fs, generally attributed to the thermalization of electrons with SCPs[16j, and a slower decay 
on the ps timescale, related to the thermalization with all other lattice vibrations. In Fig. 2C the energy-resolved 
traces at fixed delay (i=100 fs) are reported. Comparing the 5R/R{lli) measured on OP96 to the relative variation 
of the reflectivity calculated in the non-thermal and quasi-thermal cases (Fig. 1), we come to the major point of our 
work: on a timescale faster than the electron-phonon thermalization the electrons are already thermalized with some 
bosonic excitations participating to 11(17). The fast timescale (^100 fs) of this thermalization implies a very large 
coupling and a relatively small specific heat. These overall observations strongly suggest that this process involves 
bosonic excitations of electronic origin. 

The relative strengths of nb(fl) determine [2y] both the temporal evolution of the temperatures Tb, through the 
four-temperature model (4TM)[18], and the intensity of the reflectivity variation, through Eq. 1. As a consequence, 
the simultaneous fit of the calculated <5i?/i? (w, Tc,Tbo,TscPj21at) to the data reported in Fig. 2 in the time- and 
frequency domain, significantly narrows the phase-space of the parameters of the model, as compared to single-color 
measurements. The fit procedure[i8| allows us to unambiguously extract the different contributions to n(Sl) and to 
estimate Cbc and Cscp- 

Figure 3 summarizes the main results of this work. The 10 /iJ/cm^ pump pulse gently increases the electronic 
temperature by 5Tc'^2 K. The entire high-energy part and ~46% of the peak (red areas) instantaneously thermalize 
with electrons at a temperature Tbo— ^c • The spectral distribution and the value of the specific heat of these excitations 
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FIG. 3: Disentangling the contributions to the total bosonic function. The electronic (Dbc, red area), strongly-coupled phonons 
(IIscp, blue area) and lattice (Iliat, green area) contributions to the total bosonic function are reported. The insets display the 
temporal evolution of the temperatures The (red line), Tscp (blue line) and Tiat (green line), determined through the 4TM[la]. 
Sketches of the possible microscopic mechanisms at the base of the different contributions to the total I1{Q,) are reported in the 
top panels. 



(Cbe<0.1Co) demonstrate their electronic origin. On a slower timescale (100-200 fs), the electrons thermalize with 
the SCPs that represent ~20% of the phonon density of states (Cscp=0.2Ciat), but are responsible for ~34% of the 
coupling (blue area) in the peak of the bosonic function at 40-75 meV, corresponding to ~17% of the total bosonic 
function. Prominent candidates as SCPs are the buckling and breathing Cu-0 optical modes. The third and last 
measured timescale is related to the thermalization with all other lattice modes (80% of the total) that include all 
acoustic modes and the IR- and Raman-active modes involving c-axis motion of the Cu ions and provide ^20% of the 
coupling (green area) in the peak of n(ri). 

These results have important implications for the identification of the pairing mechanism in cuprates. The 
electron-boson coupling Ab=2/nb {^) /^ dil is calculated for each subset b of the bosonic excitations, considerin 
the experimental uncertainties. In the strong-coupling regime (Ab<l-5), the critical temperature can be estimated[l 
through an extended version of the Mc Millan's equation[23], containing the logarithmic-averaged frequency (f2) 
and the portion of nb(r2) that contributes to the d-wave pairing[2^. The maximal critical temperature attainable 
is calculated assuming that each nb(r2) entirely contributes to the d-wave pairing. The coupling with SCPs 
('^SCP=0-4±0.2) is in complete agreement with the values measured on similar materials via different techniques, 
such as time- resolved photoemission spectroscopy [ig, time-resolved electron diffraction flT'l and single-color high- 
resolution time- resolved reflectivity [231 ■ Although this value is rather close to the threshold of the strong-coupling 
regime|30ll3l|. the small value of fl gives rc=2-30 K, that is far from being able to account for the high-temperature 
superconductivity of the system. The coupling of electrons with all other lattice vibrations is even smaller in 
strength (Aiat=0.2±0.2) and provides an upper bound of the critical temperature of T'c~12 K. Finally, the large 
coupling constant (Abe=l-1±0.2) and the larger Cl value of the electronic excitations give Tc=105-135 K and hence 
accounts alone for the high-critical temperature. We remark that, while nscp(^) and niat(ri) are expected to be 
temperature- independent, nbs(^) increases as new magnetic excitations emerge when approaching Tc, particularly 
in the pseudogap phase [10|. Therefore, the use of nbc(i^) determined at T=300 K to estimate the QP-boson 
coupling and Tc underestimates the electronic contribution to the pairing, further supporting our conclusion of a 
dominant electronic mechanism in the superconductivity of cuprates. All the Ab values, maximum attainable criti- 
cal temperatures, and important parameters for each subset of the total bosonic function are reported in the Table S2. 



The measured values of Abe, Ascp and Aiat and the spectral distribution of the bosonic excitations strongly 
indicate that the antiferromagnetic spin fluctuations 0, Q and the loop currents Q are the most probable mediators 
for the formation of Cooper pairs. An isotope effect in the dispersion of nodal QPs has been recently reported on 
optimally-doped Bi2212, by ARPES measurements [32]. These data have been carefully analvzed|33|. demonstrating 
that the nodal isotope effect can be explained assuming that the QP-phonon coupling represents about 10% of the 
total contribution of other bosonic excitations. From Figure 3, we estimate that the contribution of SCP is less 
than 20% of the total bosonic function. This demonstrates that our results fully explain the observed isotope effect, 
eventually providing a consistent interpretation of the most important experimental results about the QP-boson 
coupling in cuprates. 

We remark that our conclusions are rather independent of the assumption of the histogram-like form of n(J7) and 
robust against modifications of the details of the equilibrium dielectric function. In fact, the outcome of this work 
strongly supports the factorization of the self-energy at T=300 K into a temperature-dependent kernel function and 
the glue function 11(51) (see Eq. 1), even under non-equilibrium conditions. Although we do not exclude a-priori 
that the upper limit used in the determination of the bosonic function can hide possible contributions to n(J7) even 
above 1 eV and that the electron-phonon coupling may cooperate in driving the superconducting phase transition, 
we demonstrate that bosonic excitations of electronic origin are the most important factor in the formation of the 
superconducting state at high temperatures in the cuprates. 

The present findings pave the way for the investigation of the electron-boson coupling in a variety of complex 
materials, ranging from transition-metal oxides to iron-based superconductors. 
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Materials and methods 

Time-resolved spectroscopy 

The dynamics of the non-equilibrium optical response is probed by ultrashort supercontinuum light pulses obtained 
by focusing the output beam of a cavity-dumped Ti:sapphire oscillator in a photonic crystal fiber (25). The spectral 
window explored by the white pulses (1.2-2 eV) is extended in the infrared region down to 0.5 eV, using an Optical 
Parametric Amplifier (OPA), seeded by a 250 kHz amplified Ti:sapphire oscillator. The pump fluence is always 
set to 10±3 /iJ/cm^. The temporal resolution of the probe pulses varies with the wavelength from 180 fs at 800 
nm to <100 fs at 2000 nm. The diameters of the pump and probe beams on the sample are 38±2 and 20±3 /zm 
(Full- Width- at-Half-Maximum) . 

Critical temperature 

In the strong-coupling formalism, the critical temperature for d-wave pairing in a Fermi liquid with pseudopotential 
li*—0 is approximately given by {27, 28): 

Tc == 0.83iiexp[-1.04(l + Xb)/gXb] (SI) 

where \nCl—2/Xi,L Ili,{il)\nil/ndn, Xb=2 J Hi, (il) / il dil is the electron-boson coupling constant and ge[0,l] is a 
parameter that accounts for the d-wave nature of the superconducting gap. The upper bound g=l corresponds to the 
case in which Ili,{fl) entirely contributes to the d-wave pairing. 

TABLE L The important parameters that have been determined in this work for the pairing mediated by the different kinds 
of bosonic excitations are reported. The error bars of the values have been evaluated accounting for all the uncertainties in the 
experimental parameters. 
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Tc(K) 


SCPs 0.4±0.2 60 




0.2Clat 
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Lattice 0.2±0.2 47 




O.SClat 


0-12 


Electronic 1.1±0.2 87 




<0.1Ce 


105-135 



Total 1.7 69 137 



Four temperatures model 

A four-temperatures model can be used to represent the following physical processes: a short laser pulse, with 
power density (absorbed) p, impulsively increases the effective electronic temperature of the electrons with a specific 
heat Ce=^eTe (7e=7r^A^c-^(e_F)feb/3), Nc being the number of cells in the sample and N{Ep) the density of states 
of both spins per unit cell). Te will then relax through the energy exchange with all the coupled degrees of freedom 
that linearly contribute to the total Il{il)=Ili,e{^)+Ilscp{^)+^iat{^)- The rate of the energy exchange among the 
different populations is given by (26): 

dTf. _ G{Ili,e,The,Te) G(nsCP ,TsCP ,Te) G(niat,Tiat,Te) p 

dt 7eTe 7eTe 7eTe JeTe 

dTbe GiUbe,Tbe,Te) 



dt Cbe 

dTscp _ G{TlscPi TscPi Te) 
dt Cscp 

dTiat _ G{Illat,Tlat,Te) 



(53) 



(54) 



(55) 



where 



G{iib,Th,Te) = 4z2 / dnnb{n)n^[N{n,Tb) - N{n,Te)] (56) 

with b=be,SCP,lat and N{Q.,T)^{e^/*'^^ — 1)~^ the Bose-Einstein distribution at the temperature T. The specific 
heat {Cscp) of SCPs is proportional to their density of states and is taken as a fraction / of the total specific heat, 
i.e., Cscp—fCiat- We pinpoint that the 4TM used does not contain any term for the anharmonic decay of SCP {16), 
that directly couples Eqs. S4 and S5. Including this term, which is of difficult evaluation, would eventually lead to a 
smaller value of Iiiat{^) and a larger value of Iibe{^)- 

The four temperature model relies on the assumption that quasi-thermal electronic and bosonic distributions, 
thermodynamically defined by the effective temperatures T^ and Tf, (with b=be,SCP,lat), are rapidly established. 
This corresponds to assume that the electron-electron scattering time {Tg-e) is much smaller than the electron-phonon 
scattering time {Tg-ph)- The evaluation of T^-e in strongly-correlated systems and in cuprates, in particular, is a 
subtle and open problem, since the applicability of arguments based on the Fermi-liquid theory is questionable {34)- 
In Refs. {31, 35) an analytical approach, based on the Boltzmann equation and free from any quasi-equilibrium 
approximation, has been used to extract the electron-phonon coupling from single-color time-resolved measurements, 
obtaining results fully compatible with our conclusions (A^O.5 in Lai.85Sro.i5Cu04 and ^0.25 in YBa2Cu306.5). 
Using the solution of the Boltzmann equation, the authors of {31) demonstrated that the quasi-thermal models 
overestimate Xe-ph of a factor 8/5, in the extreme case Te_e/Te-p/i— >oo. In our case, the fit of the 4 temperature 
model to the time-resolved data gives Ascp=0.37 that, divided by 8/5, becomes 0.23, well within the error bars 
(Ascp=0.4±0.2) reported in our work. In conclusion, the large error bars associated to the values of Xscp and Xiat 
include: i) the experimental uncertainty in the pump fluence; ii) the possibility of adding to the 4TM a term that 
accounts for the anharmonic coupling of SGP and the lattice; iii) the possibility of the failure of the 4TM in the 
extreme case Te-e^'Te-ph- 



Equilibrium dielectric function 

Figure SI reports the afe-plane reflectivity, as measured at 300 K on optimally doped Bi2Sr2Cao.92Yo.o8Cu208+5 
(Y-Bi2212) high-quality crystals (Tc=96 K), by conventional spectroscopic ellipsometry {21). The dielectric function 
has been obtained by applying the Kramers-Kronig relations to the reflectivity for 50 cm~^<ci;/27rc<6000 cm~^ and 
directly from ellipsometry for 1500 cm~^<a;/27rc<36000 cm~^. This combination allows a very accurate determination 
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FIG. 4: The real part (ei, right axis) and the imaginary part (e2, left axis) of the equilibrium dielectric function of an optimally 
doped Y-Bi2212 crystal, measured over a wide energy range at T=300 K, are reported. The inset shows the bosonic function 
n(r2), extracted from the optical spectra. The blue arrow indicates the dressed plasma frequency f2p ~ 1 eV. The black curves 
are the best fit of the extended Drude model to the data. 

of e(aj) in the entire combined frequency range. Due to the ofF-normal angle of incidence used with ellipsonietry, the 
a6-plane pseudo-dielectric function had to be corrected for the c-axis admixture. 

The solid black line in the Fig. SI is the fit to the data of an extended Drude model, describing the strong coupling 
of the free carriers with the bosonic excitations, and a sum of Lorentz oscillators accounting for the interband optical 
transitions in the visible region (19). 

Within the extended Drude model, a simplified expression of (j{uj) is obtained (32) in terms of the one-particle 
self-energy, by omitting the vertex corrections (Migdal approximation): 



aiuj,T) 



fiLU + e,T)-fie,T) 
i^TTLu Jo w - S (cj + e, T) + S* (e, T) 



de 



(57) 



where ujp is the plasma frequency and / (e, T) the Fermi distribution at temperature T. E (w) and S* (w) are the 
electron and hole self-energies (obtained by averaging S (w, k) over the Fermi surface and assuming a constant density 
of states) {24): 
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E(w,T) = / U{n)L{uj,n,T) 
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dn 



(SS) 



where L {lu, 51, T) is a material-independent kernel function accounting for the thermal excitations of the glue and the 
QPs. 

The normal-incidence reflectivity R{T, to) is calculated from Eq. S7 through the relation: 



RiT,Lu) 



i-vM7>) 



l+^e{T,iu) 



{S9) 



IKJI) is directly determined fitting the reflectivity obtained from Eq. S7-S9 to the optical spectra reported in the Fig. 
SI. 



Differential fitting procedure 



The relative reflectivity variation, i.e. SR/R{uj,t), is a functional oiIlb{^) {b=be,SCPJat) (see Eqs. S7, S9) and of 
Cf,, that are the parameters determining the temporal evolution of Tf, through the four-temperature model described 
in the previous section of the Methods. The different contributions to 11(51) are extracted fitting SR/R{uj,t) to the 
time- and frequency-resolved refiectivity data, reported in Figure 2. 



Considering that the energy distribution of phonons is hmited to <90 meV, we assume that, for f2>90 meV, 
n(ri)~nte(f^)- within this assumption, the functional dependence oiSR/R{ui, t) on Ilb{^) is simphfied as a parametric 
dependence on the coefficients pb, where \lbe{^)=Pbe^{^<^^ mcV)+n(ri>90 meV), Iiscp{^)=Pscp^{^<^^ meV) 
and I].iat{^)^Piat^{^<QQ meV). 

Considering the constraints given by the relations SR/ R{lu, t)=F{pbe^pscPiPiat, Cbe, CscPi Ciat] ^, t) (F being the 
generic function expressing the dependence on the parameters contained in Eqs. S7-S9) and n(ri)=^{,nh(il) (n(f7) 
being the total glue function reported in the Figure 2), the free parameters of the fitting are four. From the fit to 
the rise time of the time-resolved signal, it is possible to obtain an upper limit to the value of the specific heat of the 
electronic excitations, i.e., Cbe<0.1Ce. Fixing the values Ce/Te^je—W^'^ Jcm^'^K^^ and C(ai=2.27 Jcm^'^K^^, we 
are able to simultaneously perform the fit to the data reported in the Figure SI in the time and frequency domain 
with only two free parameters, i.e. Cscp &iid any of the p^- Furthermore, the small infiuence of the variations of 
Cscp on the fit results (i.e. dF/dCscP'^dF/dpb), significantly narrows the parameters phase-space of the model, 
allowing to unambiguously haul out the different contributions to n(il). 

To carry out a quantitative and reliable analysis of the results, both the finite lateral dimension of the probe beam 
and the inhomogeneous excitation of the pump pulse, related to the finite penetration depth, have been accounted for. 
Considering a pump pulse of energy P and a probe pulse with gaussian shapes in the radial coordinate, the measured 
reflectivity variation is given by: 

SR f°" 1 ^2 ,2 P ^2 ,2 P , ^ 

oc / -^e-"^ /"p-^e-" ^"^-rdr = — ^ — (SIO) 



<^pr(pu) =FWHM/2 Vln2 being the lateral size of the probe(pump) beam. As compared to the reflectivity variation 

calculated for a homogeneous excitation with a surface density of 4P/FWHM^, the measured 5R/R{^) is suppressed 
of a factor ln2-FWHMp„V(FWHMp„2+FWHMp^2) 

The role of the finite penetration depth of the pump pulse (dpu^lGO nm @ 1.55 eV) is accounted for numerically 
calculating 5R/ R{yi) through a transfer matrix method, when a graded index of the variation of the refractive index 
n with exponential profile, i.e. (5n=(5noe~^'''''" {dpu being the penetration depth of the pump pulse), is assumed. 
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